Dipolar Interactions and Origin of Spin Ice in Ising Pyrochlore Magnets 
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Recent experiments suggest that the Ising pyrochlore magnets Ho2Ti207 and Dy2Ti2C>7 display 
qualitative properties of the spin ice model proposed by Harris et al. Phys. Rev. Lett. 79, 2554 
(1997). We discuss the dipolar energy scale present in both these materials and consider how 
they can display spin ice behavior despite the presence of long range interactions. Specifically, we 
present numerical simulations and a mean field analysis of pyrochlore Ising systems in the presence 
of nearest neighbor exchange and long range dipolar interactions. We find that two possible phases 
can occur, a long range ordered antiferromagnetic one and the other dominated by spin ice features. 
Our quantitative theory is in very good agreement with experimental data on both Ho2Ti2 07 and 
Dy2Ti207. We suggest that the nearest neighbor exchange in Dy2Ti207 is antiferromagnetic and 
that spin ice behavior is induced by long range dipolar interactions. 



An exciting development has occurred in the last two 
years with the discovery of an apparent analogy between 
the low temperature physics of the geometrically frus- 
trated Ising pyrochlore compounds Ho2Ti207 jl| and 
Dy2Ti2C>7 (so called 'spin ice' materials), and pro- 
ton ordering in real ice The magnetic cations Ho 3+ 
and Dy 3+ of these particular materials reside on the py- 
rochlore lattice of corner sharing tetrahedra. Single-ion 
effects conspire to make their magnetic moments almost 
ideally Ising-like, but with their own set of local axes. In 
particular, each moment has its local Ising axis along the 
line connecting its site to the middle of a tetrahedron to 
which it belongs (see inset of Fig. |l|). 

In a simple model of nearest neighbor ferromagnetic 
(FM) interactions, such a system has the same 'ice rules' 
for the construction of its ground state as those for the 
ground state of real ice ||[|]. In both cases, these rules 
predict a macroscopically degenerate ground state, a fea- 
ture that a number of geometrically frustrated systems 
possess ^-^]. 

In Ho2Ti207, /iSR data indicates a lack of ordering 
down to ~ 50 mK despite a Curie- Weiss temperature 
#cw ~ 1-9 K, while single crystal neutron scattering data 
suggests the development of short-range FM correlations, 
but the absence of ordering down to at least 0.35 K |Q. 
Ho2Ti207 also displays field dependent behavior consis- 
tent with a spin ice picture j^] . Quite dramatically, ther- 
modynamic measurements on Dy2Ti20y show a lack 
of any ordering feature in the specific heat data, with 
the measured ground state entropy within 5% of Paul- 
ing's prediction for the entropy of ice S . 

However, both spin ice materials contain further in- 
teractions additional to the nearest neighbor exchange. 
Often, rare earth cations can have appreciable magnetic 
moments and, consequently, magnetic dipole-dipole in- 
teractions of the same order as, if not larger than the 
exchange coupling, can occur. Furthermore, it has been 
suggested that the nearest neighbor exchange interac- 
tion in Ho2Ti207 is actually antiferromagnetic (AF) [fl0|| , 



which by itself should cause a phase transition to a long 
range ordered ground state. Thus, how these systems ac- 
tually display spin ice-like behavior is most puzzling. For 
example, one might naively expect that the long-range 
and anisotropic spin-space nature of the dipolar interac- 
tions would introduce so many constraints that a large 
degree of the degeneracy present in the simple nearest- 
neighbor ferromagnetic spin-ice model |l],|j would be re- 
moved, and induce long-range order. It is this issue that 
we wish to address in this paper. 
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FIG. 1. Phase diagram of Ising pyrochlore magnets with 
nearest neighbor exchange and long range dipolar interac- 
tions. Dnn and J nn are the Ising parameters for the near- 
est neighbor dipolar and exchange energies respectively (see 
text). Inset: An ice rule configuration of a tetrahedron. The 
four local (111) Ising axes meet in the middle of the tetrahe- 
dron [4]. 



A previous attempt to consider dipolar effects in Ising 
pyrochlores was made by Siddharthan et al. |l(| . In that 
work, the dipole-dipole interaction was truncated beyond 
five nearest neighbor distances, and a sharp transition 
between paramagnetism and a partially ordered phase 
(where rapid freezing occurs) was observed for interac- 
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tion parameters believed appropriate for Ho2Ti207 [ [10| . 
However, we argue that the truncation of dipole-dipole 
interactions can be misleading, and may introduce spu- 
rious features in various thermodynamic properties. For 
example, we find that the sharp feature observed in the 
specific heat for truncation beyond five nearest neighbor 
distances JTo| ] is softened and rounded for truncation be- 
yond the tenth nearest neighbor shell, and the observed 
dynamical freezing is pushed to lower temperatures. As 
we show below, in the limit of infinite range dipoles, the 
interaction parameters of Siddharthan et al. JTcfl yield 
spin ice. 

In this work, we consider the interplay between near- 
est neighbor exchange and dipolar interactions by taking 
into account the long range (out to infinity) nature of 
the dipolar interactions through the use of Ewald sum- 
mation techniques fll|]i2"|| . Our Monte Carlo simulations 
and mean field results show that dipolar forces are re- 
markably adept at producing spin ice physics over a large 
region of parameter space. 

Some of our main conclusions are shown in Fig. [l]. 
For Ising pyrochlores, the dipole-dipole interaction at 
nearest neighbor is FM, and therefore favors frustration. 
Beyond nearest neighbor, the dipole-dipole interactions 
can be either FM or AF, or interestingly, even multiply 
valued [l3| , depending on the neighbor distance. Defin- 
ing the nearest neighbor dipole-dipole interaction as D nn 
and the nearest neighbor exchange as J nn , our Monte 
Carlo results indicate that spin ice behavior persists in 
the presence of AF exchange up to J nn /D nn ~ —0.91. 
For J, m /D nn < —0.91 we find a second order phase tran- 
sition to the globally doubly degenerate Q = phase of 
the nearest neighbor AF exchange-only model jl5| , where 
all spins either point in, or all out of a given tetrahedron. 

Our Hamiltonian describing the Ising pyrochlore mag- 
nets is as follows, 
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where the spin vector S^* labels the Ising moment of mag- 
nitude \S\ — 1 at lattice site i and local Ising axis Zi. 
Because the local Ising axes belong to the set of (111) 
vectors, the nearest neighbor exchange energy between 
two spins i and j is J nn = J/3. The dipole-dipole in- 
teraction at nearest neighbor is D nn = 5D/3 where D 
is the usual estimate of the dipole energy scale, D = 
(^ /47r) 5 2 /i 2 /^nn- For both Ho 2 Ti 2 7 and Dy 2 Ti 2 7 , 
D nn ~ 2.35 K. 

It is well known in the field of electrostatic interactions 
that the dipole-dipole interaction is difficult to handle 
due to its 1/r 3 nature. In general, a lattice summation 
of dipole-dipole interactions is conditionally convergent, 



and must be considered with care. In order to include the 
important long range nature of the dipole-dipole interac- 
tion, we have implemented the well known Ewald method 
within our simulation technique, in order to derive an ef- 
fective [ |l4| dipole-dipole interaction between spins within 
our simulation cell [ ]l2] j. Unlike in dipolar fluid simula- 
tions, our lattice spins allow the effective interactions be- 
tween all moments to be calculated only once, after which 
a numerical simulation can proceed as normal. 

A standard Metropolis algorithm was used in our 
Monte Carlo simulations. We used a conventional cubic 
cell for the pyrochlore lattice, which contains 16 spins. In 
general, we found it sufficient to simulate up to 4 x 4 x 4 
cubic cells (ie. L=4, and 1024 spins) with up to ~ 10 6 
Monte Carlo steps per spin when necessary. Thermody- 
namic data were collected by starting the simulations at 
high temperatures and cooling very slowly. 

Referring to Fig. El the characterization of a system as 
having spin ice behavior was carried out by determining 
the entropy, via numerical integration of the specific heat 
divided by temperature data. Pauling's argument for the 
entropy of ice yields R[ln 2- (1/2) ln(3/2)] or 4.07 J moh 1 
K- 1 §,§. We find for J„ n /Am = -0.91, (our spin ice 
data point closest to the phase boundary in Fig. pj) , this 
value for the entropy to within 3%, using a system size 
L = 4. 

Our thermodynamic data indicates that when the 
nearest neighbor exchange is AF and large compared to 
the dipolar interactions, the system undergoes a second 
order phase transition to an all in or all out Q — ground 
state as alluded to earlier. This AF phase persists slightly 
beyond the point where the nearest neighbor dipolar in- 
teraction (FM) is stronger than nearest neighbor AF ex- 
change. 
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FIG. 2. Dependence of the specific heat peak height and 
temperature location on exchange and dipole-dipole interac- 
tion parameters. 



In the spin ice regime, our specific heat data has a num- 
ber of interesting features which help shed light on the 
effect of long range dipole-dipole interactions. In general, 
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each data set shows qualitatively the same broad peak 
as observed in the nearest neighbor FM exchange model, 
and vanishes at high and low temperatures ||[l(| . We de- 
tect very little system size dependence upon comparison 
of data from L = 2, 3,4 simulation sizes. In particular, 
there is no noticeable size dependence of the specific heat 
maximum, nor its position. As the AF/spin ice phase 
boundary is approached from the spin ice side, the spe- 
cific heat peak begins to narrow and more importantly, 
both the peak height and peak position begin to vary. 
As we discuss below, this has important ramifications for 
the interpretation of experimental data. 

In Fig. H we plot the dependence of the specific heat 
peak height (C pea k) and the temperature at which it oc- 
curs (Tp 0a k) on the ratio of the nearest neighbor exchange 
and dipolc-dipole energies, Jnn/Dnn- Note that in the 
regime of large nearest neighbor FM coupling, the peak 
height plateaus to the value one observes in the nearest 
neighbor FM exchange-only model. 

Indeed our data suggests that when the exchange be- 
comes FM, the nearest neighbor effective bond energy is 
large enough to dominate the excitations of the system. 
This can be more dramatically seen in Fig. ^ where we 
have rescaled the specific heat curves for a number of in- 
teraction parameter values in terms of the effective near- 
est neighbor interaction J G ff = J nn + D nn in the regime 

Jnn/Am > 0. 

3.0 -i 1 1 1 1 1 1 1 1 1 1 
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FIG. 3. Specific heat for system size L = 2 rescaled 
into units of the effective nearest neighbor interaction J e g. 
Jnn/Dnn = corresponds to purely dipolar interactions while 
Jun/D nn = oo corresponds to nearest neighbor FM exchange 
only. 



short range truncation of the dipole-dipole interaction. 

As the nearest neighbor exchange interaction becomes 
AF, we find that the approximate 'collapse' onto a single 
energy scale becomes less accurate, with features in the 
specific heat becoming dependent on J nn /D nn in a more 
complicated manner. It is within this regime that we 
believe that both Ho2Ti207 and Dy2Ti2C>7 exist. 

Since D nn is a quantity which can be calculated once 
the crystal field structure of the magnetic ion is known, 
the nearest neighbor exchange J nn is the only adjustable 
parameter in our theory. Fig. ^| enables us to test in two 
independent ways the usefulness of our approach to the 
long range dipole problem in spin ice materials. 

If we consider Dy2Ti2C>7 for example, specific heat 
measurements by Ramirez et al. [0 indicate a peak 
height, Cp e y ak of - 2.72 J mol^R- 1 . Given that 
D nn « 2.35 K for this material, the left hand plot of 
Fig. ^| indicates a nearest neighbor exchange coupling 
Jnn ~ —1.2 K. The same experimental specific heat 
data shows that this peak occurs at a temperature of 
T pcak D Y ~ 1.25 K. Using the plot of T poa k m Fig. ||, wc 
independently arrive at approximately the same conclu- 
sion for the value of the nearest neighbor exchange. Thus, 
we predict that AF exchange is present in Dy 2 Ti 2 07 with 
Jnn ~ —1-2 K. If there was no AF exchange present in 
this system, our results in Fig. [j] imply that there would 
be a peak in the specific heat at a temperature of at least 
~ 2.3 K, which is not observed experimentally ||. Our 
best fit for the specific heat data of Dy2Ti2C>7 by Ramirez 
et al. Q is shown in Fig. ||, where we find good agree- 
ment between theory and experiment for J nn = —1.24 
K. 



This figure shows that in terms of an effective energy 
scale, the medium to long range effects of the dipolar in- 
teractions are 'screened' by the system, and one recovers 
qualitatively the short range physics of the nearest neigh- 
bor spin ice model. Remarkably, inclusion of long range 
dipolar interactions appears to have the effect of remov- 
ing a tendency towards ordering which can come from a 
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FIG. 4. Comparison of (a) specific heat and (b) entropy 
data between Dy2Ti207 [2] and Monte Carlo simulation with 
J nll = —1.24 K, _D nn = 2.35 K and system size L=4. 

Specific heat measurements on a powdered sample of 
Dy2Ti207 in a magnetic field were also reported in Ref. 
0. Three field independent peaks were observed at low 
temperature. For a large field in the (110) direction, two 
spins on each tetrahedron are pinned by the field, while 
the other two remain free, since their Ising axes are per- 
pendicular to the applied field. Due to the dipolar inter- 
action, there will be a coupling between the fluctuating 
spins on these two sub-lattices. Our preliminary simu- 
lations on small lattice sizes suggest a field independent 
ordering at low temperature as observed in experiment. 

Considering Ho2Ti207, experimental data on its ther- 
modynamic properties is not so categorical. The specific 
heat data of Siddharthan et al. jl0| indicates a feature at 
~ 0.8 K, although it has been suggested that this could 
be due to an additional contribution in this temperature 
range from an anomalously large hyperfine coupling in 
Ho 3+ HJ. Nevertheless, using the plot of T poa k in Fig. 
||, we find substantial AF exchange coupling of the same 
order of magnitude as in the Dy compound. We note that 
Siddharthan et al. pQ] and den Hertog et al. |D| find from 
analysis of magnetization measurements a similar order 
of magnitude for J nn . Furthermore, our numerical simu- 
lations within this region of parameter space indicate a 
Curie- Weiss temperature of ~ 2 K, in agreement with an 
experimental estimate by Harris et al. of 9 CW ~ 1.9 K [Q. 

Tb2Ti207 is an Ising pyrochlore system of similar type 
to the Ho 3+ and Dy 3+ based materials [0J|], but the 
Ising anisotropy is reduced to much lower temperature 



due to narrowly spaced crystal field levels || . While this 
makes the interpretation of experimental data more dif- 
ficult, initial estimates of the nearest neighbor exchange 
and dipole moment Q yield J n n/D nn ~ — 1, placing this 
system very close to the phase boundary of Fig. [l| In- 
deed, fiSR measurements suggest that Tb2Ti 2 07 fails to 
order down to 70 mK and thus it would appear that a 
spin ice picture for this material cannot be a priori ruled 
out. 

While we believe our approach yields a reasonably suc- 
cessful quantitative theory of spin ice behavior in Ising 
pyrochlores, there still remains the question of why long 
range dipolar interactions do not appear to lift the macro- 
scopic degeneracy associated with the ice rules, and select 
an ordered state. 

Mean-field theory provides a more quantitative basis 
for examining this issue. Following the approach used 
in Refs. [ ^6|JT7[| , we find that the soft-modes for Ising 
pyrochlore systems described by Eq. [j] consist of two 
very weakly dispersive branches (less than 1% disper- 
sion) over the whole Brillouin zone (except at Q = 0). 
Such a set of quasi-dispersionless branches is very similar 
to the two completely dispersionless soft branches of the 
nearest-neighbor FM spin ice model. Consequently, both 
nearest-neighbor FM and dipolar spin ice behave almost 
identically over the whole temperature range spanning 
0(w) < T < oo where w is the bandwidth of the soft 
branch of dipolar spin ice. We note that this near disper- 
sionless behavior is only recovered asymptotically as the 
long range dipoles are included out to infinity. The lifting 
of degeneracy at the 1% level in the apparent absence of 
any small parameters in the theory is at this point not 
completely understood. We do not believe it is due to 
numerical or computational error. A partial explanation 
may be that (111) Ising anisotropy, the 1/r 3 long-range 
nature of dipolar interactions, and the specific relation- 
ship between the topology of the pyrochlore lattice and 
the anisotropic (spin— space) coupling of dipolar interac- 
tions, combine in a subtle manner to produce a spectrum 
of soft modes that approximate very closely the spectrum 
of the nearest-neighbor ferromagnetic spin ice model of 
Harris et al. In other words, there must be an almost 
exact symmetry fullfilled in this system when long-range 
dipolar interactions are taken into account. However, 
the same long-range nature of these interactions renders 
it difficult to construct a simple and intuitive picture of 
their effects, and we have not been able to identify such 
"almost exact" symmetry. 

Also, in these systems, the absence of soft fluctuations 
(Ising spins) combined with the macroscopic degeneracy 
associated with the ice rules may be such that correla- 
tions associated with a 'true' ground state are dynam- 
ically inhibited from developing. For energetic reasons, 
states obeying the ice rules are favored down to low tem- 
perature, by which time such large energy barriers ex- 
ist that evolution towards the true ground state is never 



4 



achieved. In simulation terms, at low temperature the 
Boltzmann weights for local spin flips 'towards' this state 
are simply too small. 

In conclusion, using Ewald summation techniques we 
have considered the effects of long range dipole-dipole in- 
teractions on the magnetic behavior of Ising pyrochlorc 
systems. Our results show that spin ice behavior is re- 
covered over a large region of parameter space, and we 
find quantitative agreement between our approach and 
experimental data on spin ice materials. 
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